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Abstract 

Promising pathways for photoassociative formation of ultracold 
heteronuclear alkali metal dimers in their lowest rovibronic levels (de- 

rn ■ 

t^J- ■ noted X(0,0)) are examined using high quality ab initio calculations 

Q^ ■ of potential energy curves currently available. A promising pathway 

■ for KRb, involving the resonant coupling of the 2 1 ]! and l 1 !! states 

just below the lowest excited asymptote (K(4s)+Rb(5pi/2)) ; is found 
to occur also for RbCs and less promisingly for KCs as well. The res- 
^ ' onant coupling of the 3 1 E + and l 1 !! states, also just below the lowest 

excited asymptote, is found to be promising for LiNa, LiK, LiRb, and 
less promising for LiCs and KCs. Direct photoassociation to the l 1 !! 
state near dissociation appears promising in the final dimers, NaK, 
NaRb, and NaCs, although detuning more than 100 cm" 1 below the 
lowest excited asymptote may be required. 



1 Introduction 



Research with ultracold atoms has been an amazingly fruitful area of frontier 
research for over two decades, with many exciting developments still occur- 
ing frequently. Ultracold research has expanded into other areas of atomic, 
molecular, and optical physics including ultracold ions, plasmas, and espe- 
cially molecules, where similarly exciting developments are now occuring [1] . 

In ultracold molecules, emphasis has shifted to ultracold polar molecules, 
e.g. heteronuclear alkali metal dimers, because of the permanent dipole mo- 
ments of these molecules, which vary from small (KRb) to large (LiCs) [2]. 
Emphasis has also shifted from study of the weakly bound levels formed 
by photoassociation (PA) (1-100 cm -1 below dissociation asymptotes) and 
extremely weakly bound levels formed by magnetoassociation (MA) via Fes- 
hbach resonances to formation of the most strongly bound levels, i.e. the 
rovibronic ground state, X X S + (v—0, J=0), or X(0,0) for short. [Actually 
this lowest "level" consists of a large number of hyperfme states with very 
small splittings not normally observable in high resolution spectroscopy [3].] 
The major motivation for seeking X(0,0) molecules is their lack of inelastic 
collisions at ultracold temperatures [although the hyperfme state distribution 
can certainly be redistributed by inelastic collisions]. A second motivation 
is that the dipole moment is expected to be largest for the X(0,0) level, an 
order of magnitude larger than for levels formed by PA and many orders of 
magnitude larger than for levels formed by MA [J]. 

The formation of X(0,0) levels in ultracold molecules was first achieved in 
K2 by spontaneous emission following two-photon PA [5] and then in RbCs 
using stimulated Raman following spontaneous emission following PA [6] and 
in LiCs by spontaneous emssion following far-detuned PA [7j. More recently, 
KRb formed by MA has been efficiently transferred by stimulated Raman to 
X(0,0) [8] and CS2 formed by PA has been optically pumped to X(0,0) [9]. 
Further exciting results are expected in the very near future. 

Two alternate approaches for forming X(0,0) molecules have recently been 
proposed: Feshbach Optimized Photoassociation (FOPA) [TU] and Resonant 
Coupling of PA Excited States [H]-[13]. It is the latter approach which 
we will explore here for all 10 heteronuclear alkali metal dimers, beginning 



with the best understood case of KRb, which we have studied extensively at 
UConn [13J. It should be noted that FOPA could be described as "Resonant 
Coupling of PA Lower States" and is closely related to our approach. 

Resonant coupling of levels in two (or more) electronic states is, of course, 
a very well known and long explored topic in electronic spectroscopy; we 
say the levels mutually perturb each other or that the levels are of mixed 
character. When levels from different electronic states have nearly the same 
energy and quantum numbers, their wavefunctions become mixtures of single 
electronic channel wavefunctions. A particularly well studied example is the 
resonant coupling of the levels of the A 1 E+ and b 3 IT u states of Na2 [14] . This 
coupling has been known since the work of R. W. Wood over a century ago. 
However, the resonant couplings we are looking for are those which: (1) occur 
in the PA energy regions (<100 cm _1 below the PA asymptote); and (2) have 
good Franck-Condon wavefunction overlap with the X(0,0) wavefunction. 
The next section describes the very promising example of KRb and then we 
shall explore the other nine heteronuclear dimers. 

2 The Example of KRb 

As discussed more completely elsewhere |13j . the potential energy curves of 
KRb are fairly well known from several sets of high-quality ab initio calcu- 
lations (e.g. Ref . [To] ) and from experiment (e.g. Refs. PS] and [17] ) . A set 
of ab initio curves including spin-orbit interactions [15] is shown in Fig{TJ 

Kato and co-workers extensively studied the l 1 !! and 2 1 n states and 
their mutual perturbations [T5].[T7]. Their studies included levels of both 
states in the PA energy region below the K(4s) + Rb(5pi/ 2 ) asymptote and 
above it; those above the asymptote were observed to predissociate. The l 1 !! 
state correlates with the higher K(4s)+ Rb(5p3/2) asymptote while the 2 1 I1 
state correlates with the much higher K(4p 3 / 2 )+ Rb(5s) asymptote. The PA 
energy region levels of both states were tabulated in Ref. [13] , and the coupled 
levels noted. Our initial assignments of our KRb PA spectra [T3] included 
observations of v=61-63 of the l 1 !! state, but these levels are not resonantly 
coupled to levels of the 2 1 n state. Recently we have identified levels v=60 
of the l 1 !! state and v=17 of the 2 l U state, which are resonantly coupled, 
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Figure 1: Potential energy curves for KRb up to the second exited atomic 
asymptote based on the data in Ref. [15J. The vertical arrows correspond 
to Photoassociation (PA) and Radiative decay. The horizontal dotted lines 
correspond to the energy of the lowest excited atomic asymptote. 

as shown in FigfSJ It is the wavefunction of the v=17 level which has good 
Franck-Condon overlap with the wavefunction of the X(0,0) level. However, 
without resonant coupling we would not have access to this level by PA. That 
is because the v=17 level is basically a short range level with an extremely 
small probability of being reached by PA, since the 2 1 n potential curve is 
dissociating to a much higher asymptote and the classically forbidden part 
of the v=17 wavefunction starts at ~11 ao (see Figl3]of Ref.[13j). However, 
since v=17 and v=60 are strongly mixed, we expect the former to provide 
good overlap with X(0,0) and the latter to provide good PA probability. 

In the following section, we examine high quality ab initio potential energy 
curves for the nine other heteronuclear alkali metal dimers [18]- [23] and dis- 
cuss their promise for similar resonant-coupling-assisted formation of X(0,0) 
dimers. This is in part due to the paucity of high quality spectroscopic results 
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Figure 2: Fragment of the PA spectrum of 39 K 85 Rb from Ref.[T3] in the 
12533.8-12535.4 cm -1 range, showing two previously assigned bands [T3] of 
the 2(1) and 3(0 + ) states (marked above) and two newly assigned bands of the 
4(1) [l 1 !!] and 5(1) [2 1 n] states (marked below), corresponding to resonantly 
coupled vibrational levels v=60 and v=17, respectively. 

for these species in the PA energy region similar to those for KRb [16], [17]. 
We also note that calculations including spin-orbit coupling are not gener- 
ally available and thus the arguments are somewhat less compelling for the 
species containing heavier alkali atoms, especially Cs. 

We note that it is more efficient to estimate strong transitions by simply 
examining the internuclear distances at which classical inner and outer turn- 
ing points occur rather than generating tables of Franck- Condon factors. 
Thus, for example, the photoassociation of CI2 was estimated by examining 
the upper state wavefunctions with good overlap with the v" wavefunctions 
[24] . as discussed in Ref . [25] ( especially in connection with Fig. 175). 

Thus to strengthen our arguments, we have plotted the Franck-Condon 
factors (FCFs) for l x n (v r ) - X 1 ^ (v"=0) transitions Fig.flSK) and 2 1 n (v f ) 
- A 1 S + (v"—0) transitions Fig.(J3]B). As can be clearly seen in FigJ3A the 
peak of the FCFs for l 1 !! (v') - A X S + (v"=0) transitions is at v'—2, which 
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Figure 3: FCFs versus internuclear distances for A) l 1 ]! (v') - (v"=0) 
transitions and B) 2 1 I1 (v') - (v"=0) transitions. . The blue coordi- 

nates are the inner classical turning points R' v _ of the upper state and the 
green are the inner turning point Rq_, the equilibrium internuclear distance 
R", and the outer turning point R'L of the X(0,0) level. All internuclear 
distances are represented on the same scale, shown at the bottom. 



would be very difficult to reach by direct PA. On the other hand, for 2 1 I1 
(V) - (v"=0) transitions, the peak is at v'=20 and v'=17 also has a 

very promising Franck-Condon overlap with X(0,0) level. However, because 
of the strong mixing (resonant coupling) of v'—17 level of the 2 1 I1 state with 
v'=60 level of the l 1 !! state and the observation of these two mixed levels in 
FigJ2J these mixed levels will each have good FCFs for emission to the X(0,0) 
level. 

3 Other Heteronuclear Dimers 

3.1 LiNajlB] 

The potential energy curve which is most promising for emission to X(0,0) 
is the 3 1 S + state (FigEJ). However, since direct PA to this short range state 
would be expected to be very weak, it is only with resonant coupling with 
the l : n lone; range state that this becomes a promising pathway. Note that 
the two states cross each other twice, at approximately 6.5 and 11 a a . 

Note that there is still an avoided crossing between the 2 x n and l*n 
states at ~7.5 a a . However, the resonantly coupled levels will be above the 
lowest excited asymptotes Li(2pj)+Na(3s). The possibility of direct PA to 
the 2 1 n state just below the Li(2s)+Na(3pj) asymptotes is not as promising 
as it might appear because the state is repulsive [T7] at large internuclear 
distance; PA through a quasibound resonance might however be possible. 

3.2 LiK [Previously unpublished calculations by Olivier 
Dulieu] 

The potential energy curve which is most promising for emission to X(0,0) 
is again the 3 X S + state (FigEJ) • Again, direct PA to this short range state 
would be expected to be extremely weak. Also resonant coupling with the 
l x n long range state makes this a promising pathway. The two states cross 
twice at ~7 and 10 a a . 

Again there is an avoided crossing of the 2 l U and l l U states, but the 
coupled levels are above the lowest excited asymptotes Li(2s)+K(4pj). 




Figure 4: Potential energy curves for LiNa up to the second exited atomic 
asymptote based on the data in Ref. [18J. The vertical arrows correspond 
to Photoassociation (PA) and Radiative decay. The horizontal dotted lines 
correspond to the energy of the lowest excited atomic asymptote. 



3.3 LiRb 19 



Again the most promising potential energy curve is the 3 1 S + short range 
state which is resonantly coupled to the l 1 !! long range state near 8.5 a Q 
(FigJU]). Note that here the effect of spin-orbit may be significant, slightly 
raising the l 1 !! state (which correlates to Li(2s)+Rb(5j9 3 / 2 ) with respect to 
the 3 1 S + state (which correlates to Li(2pi/ 2 )+Rb(5s)). 

The avoided crossing; of the 2 1 IT and l 1 !! states is no longer evident. 



3.4 LiCs|20] 

The most promising potential curve is again the 3 1 E + short range state, but 
now the resonant coupling is less clear, since the state does not cross the l 1 !! 
state in the absence of spin-orbit coupling (FigJTj). Given the large spin-orbit 
coupling in the Cs atom, it is possible that the curves will in fact cross in 
an ab initio calculation including spin-orbit. We recommend that such a 




Figure 5: Potential Energy curves for LiK up to the second exited atomic 
asymptote based on the unpublished calculations by Olivier Dulieu. The 
vertical arrows correspond to Photoassociation (PA) and Radiative decay. 
The horizontal dotted lines correspond to the energy of the lowest excited 
atomic asymptote. 



calculation be carried out, as it may help explain the very interesting results 
of Ref.m. 



3.5 NaK[21 



It appears that direct PA to the l 1 !! state is the most promising pathway in 
the case of NaK (FigJHJ). Here detunings less than 100 cm -1 may be required 
to reach X(0,0) with high probability. It does not appear that resonant 
coupling will play a significant role. 




Figure 6: Potential energy curves for LiRb up to the second exited atomic 
asymptote based on the data in Ref. [19J. The vertical arrows correspond 
to Photoassociation (PA) and Radiative decay. The horizontal dotted lines 
correspond to the energy of the lowest excited atomic asymptote. 




Figure 8: Potential energy curves for NaK up to the second exited atomic 
asymptote based on the data in Ref. [2 lj . The vertical arrows correspond 
to Photoassociation (PA) and Radiative decay. The horizontal dotted lines 
correspond to the energy of the lowest excited atomic asymptote. 




Figure 7: Potential energy curves for LiCs up to the second exited atomic 
asymptote based on the data in Ref. [20J. The vertical arrows correspond 
to Photoassociation (PA) and Radiative decay. The horizontal dotted lines 
correspond to the energy of the lowest excited atomic asymptote. 



3.6 NaRb 19 



Here again it appears that direct PA to the l l U state is the most promising, 
even more so when the l 1 !! state is shifted up by the large atomic spin-orbit 
coupling in the Rb atom (FigJH]). It does not appear that resonant coupling 
will play a significant role. 




Figure 9: Potential energy curves for NaRb up to the second exited atomic 
asymptote based on the data in Ref. [19J. The vertical arrows correspond 
to Photoassociation (PA) and Radiative decay. The horizontal dotted lines 
correspond to the energy of the lowest excited atomic asymptote. 



3.7 NaCs[20j 

Here again it appears that direct PA to the l 1 !! state is the most promising 
pathway, even more when the l 1 !! state is shifted up even more than in NaRb 
by the large atomic spin-orbit coupling in the Cs atom (Fig flQl) . It does not 
appear that resonant coupling will play a significant role. 
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Figure 10: Potential energy curves for NaCs up to the second exited atomic 
asymptote based on the data in Ref. [20J. The vertical arrows correspond 
to Photoassociation (PA) and Radiative decay. The horizontal dotted lines 
correspond to the energy of the lowest excited atomic asymptote. 



The more promising path here appears to be the 3 1 S + ~ l 1 !! resonant cou- 
pling as in the LiM dimers (3.1-3.4) (Fig JTTT) . However, the overlap is not 
as good and the reduced mass is much larger, so the extent of the X(0,0) 
wavefunction is much smaller, further diminishing the Franck-Condon fac- 
tors. Since the spin-orbit interaction is much larger in Cs than in K, an 
ah initio calculation with spin-orbit included may reduce the 2 1 IT - l 1 !! gap 
and make that resonant coupling pathway rather promising as well. Thus a 
calculation involving spin-orbit coupling is desirable here. 



3.8 KCs 22 




Figure 11: Potential energy curves for KCs up to the second exited atomic 
asymptote based on the data in Ref. [22J. The vertical arrows correspond 
to Photoassociation (PA) and Radiative decay. The horizontal dotted lines 
correspond to the energy of the lowest excited atomic asymptote. 



3.9 RbCs[23j 

As in KRb, 2 1 I1 resonant coupling with l 1 !! appears to be quite a promising 
path (Fig JT2"]) . However, a calculation involving spin-orbit coupling would 
also be particularly desirable in this case. 




Figure 12: Potential energy curves for RbCs up to the second exited atomic 
asymptote based on the data in Ref. [23J. The vertical arrows correspond 
to Photoassociation (PA) and Radiative decay. The horizontal dotted lines 
correspond to the energy of the lowest excited atomic asymptote. 



4 Summary and Discussion 

The results of the examination of FigHJ{T2] are summarized in TableHJ Our 
initial example of resonant coupling of the 2 1 I1 and l 1 !! states in KRb also 
appears relevant to the case of RbCs and possibly KCs. The resonant cou- 
pling of the 3 1 S + and l 1 !! states appears to be more relevant in the cases 
LiNa, LiK, LiRb, and possibly LiCs and KCs. Simple direct PA appears most 
promising in NaK, NaRb and NaCs, without any need for resonant coupling. 

For the direct PA cases (NaK, NaRb and NaCs), experiments as a function 
of the detuning of the PA laser to the l 1 !! state while detecting the X(0,0) 
level should readily identify the optimum l 1 !! vibrational level for forming 
X(0,0) molecules. 

For the 2 1 n ~ l 1 !! resonant coupling cases (RbCs and KCs), spectroscopy 
comparable to that carried out for KRb remains to be carried out. Detecting 
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the X(0,0) level while tuning the PA laser should again allow ready identifi- 
cation of the resonantly coupled levels. 

For the 3 : S + ~ l 1 !! resonant coupling cases (LiNa, LiK, LiRb, and possi- 
bly LiCs and KCs), spectroscopy is again needed. Detecting the X(0,0) level 
while detuning the PA laser should again assist in identifying the resonantly 
coupled levels. 

In both cases of resonant coupling, there are two factors to consider: the 
magnitude of the coupling of the two single channel electronic states, and 
the energy denominator. The magnitudes of the short range couplings of two 
singlet states discussed in Section III probably do not vary greatly, although 
the homogeneous perturbations (AA = 0) are often somewhat larger than 
the heterogeneous ones (AA = ±1). However, we have thus far considered 
electronic states to be eigenfunctions of total electronic spin S. However, the 
long range states formed by PA are generally not eigenfunctions of S, but 
rather of J 2 , with eigenvalue Q, the projection of the total electronic angular 
momentum along the internuclear axis z. Thus for a particular lowest excited 
asymptote of a heteronuclear alkali dimer (say ns + n' P1/2), there are three 



states (0 + , CT, and 1) and for the second lowest excited asymptote (ns + 
n ' P3/2), there are five states (0 + , CT, 1, 1, and 2). While the 0~ and 2 
states are pure triplets, the other states are each singlet /triplet mixtures. 
Moreover, different states of the same Q value, e. g. the three Q=l states, 
can also resonantly couple, as can states with (AQ = ± 1). Thus the l 1 !! 
state we speak of being formed by PA (actually a 1 state correlating with the 
second lowest excited asymptote) may not be a pure singlet state for resonant 
coupling. Nevertheless, in the region inside 15 a a , it is well described as a 
singlet state. The states to which the resonant coupling occurs (2 1 I1 and 
3 1 E + ) are in fact well described as singlets as long as they are well below 
their dissociation limits (to the fourth and third lowest excited asymptotes, 
respectively), which they always are in the PA region below the lowest excited 
asymptote. 

The other factor, the energy denominator, can greatly increase the strength 
of resonant coupling when very small. The odds of very small energy denom- 
inators scale with the density of vibrational levels, i. e. with the square 
root of the reduced mass. Thus strong resonant coupling corresponding to 
very small energy denominators is most favored for RbCs and least favored 
for LiNa; however, this is just a factor of 3.3 in square root of the reduced 
mass. In 39 K 85 Rb [T3j , for example, in addition to the (60,17) pair of reso- 
nantly coupled levels shown in FigfS} there is the (64,18) pair ~2 cm -1 above 
the K(4s) Rb(5pi/ 2 ) asymptote and the (57,16) pair ~87 cm -1 below that 
asymptote. Thus the number of pairs of resonantly-coupled vibrational levels 
will probably remain in the single digits for all ten systems. 

It might also be noted that other useful comparisons among the het- 
eronuclear alkali dimers have previously been published: a comparison of the 
strength and signs of the long range potentials and the probability of PA[26J, 
a comparison of the theoretical and experimental dipole moments[2j, a com- 
parison of photoassociation and molecule formation probabilities [27] and a 
comparison of stimulated Raman schemes for forming X(0,0) molecules via 
the resonantly coupled 2 X E + ~ 1 3 I1 states starting in a magnetoassociated 
state of predominantly triplet character [28J. 
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